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We investigate the low-energy properties of antiferromagnetic quantum XXZ spin chains with couplings fol- 
lowing two-letter aperiodic sequences, by an adaptation of the Ma-Dasgupta-Hu renormalization-group method. 
For a given aperiodic sequence, we argue that in the easy-plane anisotropy regime, intermediate between the 
XX and Heisenberg limits, the general scaling form of the thermodynamic properties is essentially given by the 
exactly known XX behavior, providing a classification of the effects of aperiodicity on XXZ chains. As rep- 
resentative illustrations, we present analytical and numerical results for the low-temperature thermodynamics 
and the ground-state correlations for couplings following the Fibonacci quasiperiodic structure and a binary 
Rudin-Shapiro sequence, whose geometrical fluctuations are similar to those induced by randomness. 
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At low temperatures, the interplay between lack of 
translational invariance and quantum fluctuations in low- 
dimensional strongly correlated electron systems may induce 
novel phases with peculiar behavior In particular, random- 
ness in quantum spin chains may lead, for instance, to Grif- 
fiths phases fT], large-spin formation fz], and random-singlet 
phases |3]. On the other hand, studies on the influence of de- 
terministic but aperiodic elements on similar systems (see e.g. 
m 1^ 1^ I2]), inspired by the experimental discovery of 
quasicrystals, have revealed strong effects on dynamical and 
thermodynamic properties, but far less is known concerning 
the precise nature of the underlying ground-state phases. 

Prototypical models for those studies are spin- 5 antifeiTO- 
magnetic (AFM) XXZ chains described by the Hamiltonian 

H = Y.JiiSfSUi+StSl,+AS^Sl^), (1) 

where all 7, > and the 5, are spin operators. Random- 
bond versions of these systems have been much studied by 
a real-space renormalization-group (RG) method introduced 
iiai by Ma, Dasgupta and Hu (MDH) for the Heisenberg 
chain (A = 1) and more recently extended by Fisher (J 
who gave evidence that the method becomes asymptotically 
exact at low energies. The idea is to decimate the spin pairs 
coupled by the strongest bonds (those with the largest gaps 
between the local ground state and the first excited multiplet), 
forming singlets and inducing weak effective couplings be- 
tween neighboring spins, thereby reducing the energy scale. 
For XXZ chains in the regime — ^ < A < 1 , the method pre- 
dicts the ground state to be a random-singlet phase, consisting 
of arbitrarily distant spins forming rare, strongly correlated 
singlet pairs yi]. Here we employ the MDH method to inves- 
tigate the low-energy properties of aperiodic chains. 

Two-letter aperiodic sequences (AS's) can be generated by 
an inflation rule such as a ab, b ^ a, which produces the 
Fibonacci sequence abaababa . . .. Associating with each a a 
coupling Ja and with each b a coupling Jb we can build an ape- 
riodic quantum spin chain. In the XX limit (A = 0), the low- 



temperature thermodynamic behavior can be qualitatively de- 
termined for any AS by an exact RG method Iq]. The effects 
of aperiodicity depend on topological properties of the AS. If 
the fraction of letters a (or b) at odd positions is different from 
that at even positions (i.e., if there is average dimerization), 
then a finite gap opens between the global ground state and the 
first excited states, and the chain becomes noncritical. Other- 
wise, the scaling of the lowest gaps can be classified according 
to the wandering exponent co measuring the geometric fluctu- 
ations g related to nonoverlapping pairs of letters |6], which 
vary with the system size N as g ^ N"'. If CO < 0, aperiodicity 
has no effect on the long-distance, low-temperature proper- 
ties, and the system behaves as in the uniform case, with a 
finite susceptibility at T = 0. If CO = 0, as in the Fibonacci 
sequence, aperiodicity is marginal and may lead to nonuniver- 
sal power-law scaling behavior of thermodynamic properties. 
If CO > 0, aperiodicity is relevant in the RG sense, affecting 
the r = critical behavior and leading to exponential scal- 
ing of the lowest gaps (A) at long distances r, of the form 
A ^ exp(— r"). In particular, for sequences with (a= j, geo- 
metric fluctuations mimic those induced by randomness, and 
the scaling behavior is similar to the one characterizing the 
random-singlet phase 1 3]. No analogous results exist for gen- 
eral XXZ chains, although bosonization and density-matrix 
RG (DMRG) calculations on the Heisenberg chain indicate 
that Fibonacci couplings should be relevant |7, 8, 9]. We ar- 
gue below that, for a given AS, low-temperature properties of 
all chains in the regime < A < 1 should follow essentially 
the XX scaling form. Moreover, we also obtain information 
on ground-state correlation functions. 

In order to apply the MDH method to aperiodic chains, we 
must remember that now there are many spin blocks with the 
same gap at a given energy scale. Also, those blocks may 
consist of more than two spins, in which case effective spins 
would form upon renormalization. The strategy is to sweep 
through the lattice until all blocks with the same gap have been 
renormaUzed, leading to new effective couplings (and possi- 
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bly spins). Then we search for the next largest gap, which 
again corresponds to many blocks. When all possible original 
blocks have been considered, there remains some unrenormal- 
ized spins, possibly along with effective ones, defining new 
blocks which form a second generation of the lattice. The 
process is then iterated, leading to the renormalization of the 
spatial distribution of effective blocks along the generations; 
for small enough coupling ratios (which are indeed required 
for the MDH method to work), this distribution will be the 
same for allQ <A< 1 . Because of the self-similarity inherent 
to AS's generated by inflation rules, the effective couplings 
take a finite number of values, and it is natural that the block 
distribution reaches a periodic attractor (usually a fixed point) 
after a few lattice sweeps. By studying recursion relations for 
the effective couplings, we can obtain analytical results. As 
the RG proceeds, the coupling ratio usually gets smaller, sug- 
gesting that the method becomes asymptotically exact. This 
picture holds for marginal (co = 0) and relevant (co > 0) aperi- 
odicity; for irrelevant AS's, such as the Thue-Morse sequence 
(a ab, b ba), the coupling ratio approaches unity as the 
RG proceeds, and the method eventually breaks down. As 
representative examples, we consider the marginal case of Fi- 
bonacci couplings and the relevant case of a binary Rudin- 
Shapiro sequence. Full details of the calculations, as well as 
application to other AS's, will be reported elsewhere fTHl . 

The blocks to be renormalized consist of n spins connected 
by equal bonds Jq, and coupled to the rest of the chain by 
weaker bonds 7/ and Jr. The ground state (GS) for blocks 
with an even number of spins is a singlet, and at low energies 
we can eliminate the whole block, along with 7/ and Jr, leav- 
ing an effective AFM bond J' coupling the two spins closer 
to the block and given by second-order perturbation theory as 
J' = jnJiJi /Jo, with A-dependent coefficients y„. A block with 
an odd number of spins has a doublet as its GS; at low ener- 
gies, it can be replaced by an effective spin connected to its 
nearest neighbors by AFM effective bonds J'l ^ = JnJi.r whose 
values are calculated by first-order perturbation theory. In 
general, the anisotropy parameters are also renormalized and 
become site-dependent; for n even, the effective anisotropy 
is A' = 5„(Ao)A/Ar, while for n odd A'l ^ — 5„(Ao)A;,r, with 
|5„(A)| < 1 forO< A< 1 and5„(l) = L Thus,forO< A< 1 
the A,- flow to the XX fixed point (all A, = 0), ultimately re- 
producing the corresponding scaling behavior, while for the 
Heisenberg chain all A,- remain equal to unity. So, we focus 
here on the Heisenberg and XX limits, and postpone examples 
for intermediate cases to a future publication 1 1 1|. 

First we apply the method to chains with Fibonacci cou- 
plings. This is the simplest example of the quasiperiodic 
precious-mean sequences with marginal fluctuations |'6]. A 
few bonds closer to the left end of the original chain, along 
with induced effective couplings, are shown in Fig. for 
Ja < Jb ifliil . Only singlets are formed by the RG process. 



FIG. 1 : Left end of the Fibonacci XXZ chain. Dashed (solid) lines 
represent Ja (Jb) bonds. An effective coupling is induced between 
spins separated by only one singlet pair, while J'^ connects spins sep- 
arated by two singlet pairs. Apart from a few bonds close to the chain 
ends, the effective couplings also form a Fibonacci sequence. 



producing two different effective couplings, 

J'a^iiJl/jl and J'k^yiJl/Jb 



The bare coupling ratio is p = Ja/Jb^ its renormalized value 
being p' = Y2P- In each generation j, all decimated blocks 
have the same size rj and gap Aj (proportional to the effective 
Jb bonds). The recursion relations for p and A are given by 

Pj+i=y2Pj and A;+i =Y2p]A;. 

The distance between spins forming a singlet in the jth gen- 
eration defines a characteristic length rj, corresponding to the 
Fibonacci numbers rj — 1, 3, 13, 55, . . .; for j 1 the ratio 
rj+i/rj approaches (|)^, where (|) = (1 + \/5)/2 is the golden 
mean. So we have rj ^ ro(|)^', where ro is a constant, and by 
solving the recursion relations we obtain the dynamic scaling 
behavior. 
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with ^ = — |lnp/ln(|) and jj ~ — lny2/91n^(|). For the Heisen- 
berg chain Y2 = 5, and Eq. Q describes a weakly exponential 
scaling, but not of the form A ~ exp(— r") found for the XX 
chain with relevant aperiodicity (co > 0) and used to fit the 
DMRG data for the Fibonacci Heisenberg chain |7, 8]. For 
the XX chain 72 = 1, so that /j = and we can identify ^ with 
a dynamical critical exponent z, whose value depends on the 
coupling ratio, leading to nonuniversal scaling behavior, char- 
acteristic of strictly marginal operators flS*]. This nonuniver- 
sality should hold in the anisotropy regime < A < 1 with a 
"bare" value of p defined at a crossover scale. Note that we 
can view the Heisenberg scaling form (p 0) as a marginally 
relevant (co 0+) case. 

The susceptibility x{T) can be estimated d by assuming 
that, at energy scale Aj ^ T, only unrenormalized spins are 
magnetically active (and essentially free), singlet pairs be- 
ing effectively frozen. Thus, if rij ^ '"j' is the number of 
surviving spins in the jth generation, x{T ^ Aj) nj+\/Aj. 
As shown in Fig. |3a), yi{T) estimated for the Fibonacci XX 
chain from the MDH method agrees very well with results 
from free-fermion fl3\ numerical diagonalization (ND) of fi- 
nite chains, even for p = |. 

As all singlets formed at the Jth generation have size rj 
and the block distribution is fixed, the average GS correlation 
between spins separated by a distance rj is 



C««(o) ^ {SfS-^\ = i |co| («, -n,+0 = a|co| rjS (3) 
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FIG. 2: Behavior of x(r) for aperiodic chains, obtained from 
the MDH method (symbols) and from numerical diagonalization on 
chains with 10'*-10^ sites (curves), for two values of p = Ja/Jh- (2) 
Fibonacci couplings. The slope of the curves depends on p, reflecting 
the marginal character of the aperiodicity. (b) Rudin-Shapiro cou- 
plings. The inset plots the inverse square root of T%{T) versus T 
(in log scale) with p = ^ , showing that for points corresponding to 
the smallest gaps in each generation the random-singlet phase result 
XiT) ~ 1 /T\r? T is reproduced (dotted line). 



where o is a constant, a = x,>',z, and cq is the GS coiTela- 
tion between the two spins in a singlet, given by cq = — | for 
the Heisenberg chain and for both (X — x and a = z in the XX 
chain. We point out that these should be the dominant cor- 
relations, and spins separated by distances other than rj are 
predicted to be only weakly correlated. This can be checked 
for the XX chain by calculating GS correlations via the free- 
fermion method. Results for chains with p < 5 (see lHHl ) 
reveal very good agreement with Eq. 0. As correlations 
in the uniform XX chain fl^ decay as C'^'^(r) ~ r^'/^ and 
C""(r) ~ r^^, dominant xx (zz) correlations in the Fibonacci 
chain are weaker (stronger) than in the uniform chain. 

Relevant aperiodicity is characterized by strong geometric 
fluctuations, and is usually induced by sequences with blocks 
having two or more neighboring strong bonds. Furthermore, 
after the first lattice sweep, more than two values of effective 
couplings may be produced. However, they all derive from the 
original pair of couplings y„ and Jy, so that, in the presence of 
a periodic attractor, it is generally possible to write recursion 
relations for an effective coupling ratio and gap having the 
forms [11] 



P7+1 



and A 



7+1 



:XP^A;, 
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where c and X are A-dependent nonuniversal constants, and £ 
(a rational number) and k (an integer) relate to the number of 
singlets involved in determining the effective couplings. We 
assume k>2,k= \ corresponding to marginal behavior, as in 
the Fibonacci case (3- If ^ characteristic length scale takes 
the form rj = roT-', with a rescaling factor T, solving the recur- 
sion relations leads to a dynamic scaling described by 



FIG. 3: Left end of the first 3 generations of the Rudin-Shapiro 
chains. Thick lines indicate strong bonds. Shaded blocks contribute 
effective spins when renormalized; white blocks form singlets. 



with ^ and /j nonuniversal constants and co = ln^/lnx._Note 
that CO has precisely the same form as the exact result 16] ob- 
tained for the XX chain with aperiodic couplings not inducing 
average dimerization. Moreover, co depends only on the topol- 
ogy of the sequence and on its self-similar properties, but not 
on the anisotropy; so, for a given AS, the scaling form in Eq. 
(|5} should be valid for any XXZ chain in the regime < A < 1 , 
with the same exponent CO. 

As an example, we consider couplings following the binary 
Rudin-Shapiro (RS) sequence, whose inflation rule for letter 
pairs is aa —>■ aaab, ab aaba, ba bbab and bb bbba. 
This generates blocks having between 2 and 5 spins. Figure|31 
shows a few bonds closer to the left end of the chain for the 
first three generations of the lattice. Blocks with more than 
3 spins are eliminated in the first sweep and do not appear in 
later generations. Both two- and three-spin blocks are present 
in the fixed-point block distribution (already reached at the 
second generation), and a hierarchy of effective spins is pro- 
duced, as depicted in Fig. |3] At the jth generation, three-spin 
blocks have size rj ~2x 4'^ ' (so that x = 4), while two-spin 
blocks have size rj/l. The first lattice sweep generates effec- 
tive couplings Ji having 8 different values, but three of them 
are enough to write recursion relations. 



-yil^JcJl/Jl. J'b 



exp(- 



' exp (- 



(5) 



and fc='Y2'Y3JaJc/Jh- 

The gap in a given generation is proportional to the effective 
/f,; defining p — fa/Jb and eliminating fc gives 

Pi+i^lip] and Aj+i='ipy^Aj, 

which correspond to the forms in Eq. (|4} with k = 2 and 
i — 1/2. So we obtain, for the whole regime < A < 1, the 
dynamical scaling form in Eq. (|5} with an exponent CO = j, 
as predicted for the XX chain, reproducing the result for the 
random- singlet phase. In Fig. |3b) we plot x{T) for the XX 
chain calculated from both ND and the MDH method, again 
obtaining very good agreement. 

For chains with RS couplings, effective-spin formation 
from three-spin blocks determines the dominant GS coiTela- 
tions. These blocks stem from both original three-spin (and 
some five-spin) blocks and unrenormalized spins. An effec- 
tive spin represents all spins in the original block via Clebsch- 
Gordan coefficients, allowing us to calculate correlations be- 
tween any two spins whose effective spins end up in the same 
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FIG. 4: GS correlations for chains with RS couplings, obtained from 
extrapolation of numerical MDH results for chains with 2'* to 2^" 
sites, (a) C"(r) (upper solid curve) and C" {r) (lower solid curve) 
for the XX chain. The dotted curve is proportional to r^^/^. (Curves 
offset for clarity.) Inset: dominant correlations, fitted by a law 
of the form rC'*-*(r) = vo +yi Inr (dashed curve), (b) C{r) for the 
Heisenberg chain (solid curve). The dotted curve is proportional to 
I/'-- 

block at some stage of the RG process. Because of the hierar- 
chical structure seen in Fig. |3] for each block renormalized at 
the jtii generation the correlation between its end spins con- 
nects a number of order 2-' original spin pairs separated by the 
same distance rj (the size of the block), yielding a contribu- 
tion gj to the average correlation in the Heisenberg chain and 
C-"{rj) in the XX chain given by a geometric series in 2y^- For 
the Heisenberg chain ly^ = ^ < 1 , and thus 

C{rj) ^ gjiij = a\cQ\rj\ (6) 

where cq is the correlation between end spins in a three-spin 
block and rij ^ a/rj is the fraction of such blocks in the jth 
generation. For the XX chain = 1, so that ^^^{rj) carries a 
logarithmic correction, 

C'\rj)^gjnj = |co| {yQ+yx\nrj)r-]\ (7) 

where yo and y\ are constants. The zz correlation between 
end spins in a three-spin block is zero, so that the dominant 
correlations correspond to spin pairs (connected through one 
of the effective end spins and the middle spin) at distances 
r'j = 4^-1 ±4^-2 ±4^-3 ± . . . ± 1, with average (r^) = A^-\ 

and are given by g'j ~ 1 /2^^ ' . We then have 

C-(r;.)^^>;■-o'|c^,|(r^)-'/^ (8) 

Equations Q and (|8} should be contrasted with the random- 
singlet isotropic result C[r) ^ r^^, indicating a clear distinc- 
tion between the ground-state phases induced by disorder and 
aperiodicity, even in the presence of similar geometric fluctu- 
ations. This is related to the inflation symmetry of the AS's, 



which is absent in the random-bond case (or in aperiodic sys- 
tems with random perturbations 1 16]). Its effects are exem- 
plified by the fractal structure of the GS correlations visible 
in Fig. 13 which displays results from numerical implementa- 
tions of the MDH method for both XX and Heisenberg chains, 
showing conformance to the scaling forms in Eqs. (|6j-(|8|l. 

In summary, we have used the Ma-Dasgupta-Hu method to 
investigate numerically and analytically the low-energy prop- 
erties of aperiodic antiferromagnetic XXZ chains. From a 
general scaling argument, we suggest that the effects of bi- 
nary aperiodicity on the whole anisotropy regime from the XX 
to the Heisenberg limits can be classified based on the same 
wandering exponent CO which is know exactly to govern the 
scaling behavior of aperiodic XX chains. We have also shown 
that ground-state correlations are dominated by characteristic 
distances related to the rescaling factor of the sequences. 

This work has been supported by the Brazilian agencies 
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Note added. While this Letter was under review, some of 
the results for the Fibonacci chain were published in a Letter 
by K. HidaLiU- 
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